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O ■ Abstract 

(N 

^ ^. We demonstrate that the R-parity breaking interactions within their current experimental upper 

U ' bounds can give rise to large mixing-induced CP asymmetry in exclusive radiative B decays that may 

^0 , be detectable in the upcoming B factories. 
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^] I Supersymmetry is widely considered to be a leading candidate of physics beyond the standard model 

. that may be realized in Nature. So much so that the search for its signals constitutes one of the main 

' physics goals of the current and future colliders. On the other hand, there is an emerging consensus that 

■ new physics might show up through detectably large CP violation in B decays Q]. In this sense, it is 

I reasonable to ask whether it can provide an indirect but an unmistakable signature of supersymmetry. 

I Since B decay modes are often masked by theoretical as well as experimental uncertainties, one way 

^ • to establish the existence of new physics is to identify some characteristic signatures of it which can never 

be reproduced within the acceptable range of the standard model (SM) parameters. CP asymmetry in 
exclusive radiative B decays offers one such sensitive probe of physics beyond the SM. More explicitly, 
(— I I consider Bq M°7, where q is either d or s quark, and M = , oj, 0, (decaying as A'*° KsTt'^) is a 

self-conjugate meson with CP eigenvalue ^ = ±1. The CP asymmetry would arise due to the interference 
between mixing and decay. If (j)^ denotes the weak phase associated with the bji sl7l (i-c Bq — > Mjl) 
decay with amplitude Al, while and Af; are corresponding quantities for ^ sr'Jr (i-c. Bq M^n), 



Oh! 



X. 

^ then the time-dependent mixing-induced CP asymmetry will be given by. 

Tit) - fit) 

acp{t) = ^ = ^Am sin(</)M - (j^L - (pn) sin(Amt), (1) 

where, Am = 2\AlAb.\/{\Al\'^ + \Aii\'^), and 4>m is the phase of Bq-Bq mixing. In the SM, this asymmetry 
is small because of the following reason. At the quark level the interference occurs between QrJr 
whose amplitude is proportional to rriq and the hermitian conjugate (through Bq-Bq mixing) of bR — >• Ql^l 
whose amplitude is proportional to mb, leading to Am oc niq/mi,. As a result, this asymmetry is ~ 1% for 
6 — > (i7 and ~ 10% for b —^ sj. The reason for this suppression is clearly the appearance of light quark 
mass niq as a pre-factor with the b^a^^qRE^'^ part (i.e. b^ — > QrHr. decay) of the effective Hamiltonian 
vis-a-vis the appearance of mf, with the bRapiu<lLF'^'''^ part (i.e. bR — > QlIr decay). Is it possible to avoid 
this suppression by going beyond the SM? In the left-right symmetric model this asymmetry can go up to 
50% 1^ and in supersymmetric model to about 80% with large sfermion mixings |^. In this paper, we try 
to answer this question in supersymmetric models without R-parity. 

In the minimal supersymmetric standard model (MSSM), gauge invariance ensures neither the con- 
servation of Icpton number (i) nor that of baryon number {B). Defining i?-parity in terms of L and B 
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as R — (— l)(3S+L+2S)^ where S is the spin of the particle, one should in a general supersymmetric model 
allow R-parity-violating (RPV) couplings Q|. i? is +1 for all SM particles and —1 for their superpartners. 
Even though any concrete evidence for the existence of RPV terms is still lacking, the recent observation of 
neutrino masses and mixings in solar and atmospheric neutrino data suggests that it would be premature 
to abandon the i-violating interactions. However, to avoid rapid proton decay one cannot simultaneously 
admit both L- and i?-violating interactions and for this reason we impose B conservation by hand. The 
L-violating superpotential can be written as (with i, j, k as generation indices) 

W^RPV = -j^^ijkLiLjEl + \'i,jf.LiQjDl + ^jLjHu, (2) 

where A^-fe = —Xjik- Here Li and Qi are lepton and quark doublet superfields, E'[ and D'l are charged 
lepton and down quark singlet superfields, and iJ„ is the Higgs superfield that gives mass to up- type quarks. 
The trilinear Ayfe-couplings and bilinear /ij mass parameters are not relevant for our purpose and from now 
on we consider only the trilinear A^jj.-induced interactions. 

Tight constraints on the sizes of these couplings have been placed from the consideration of neutrinoless 
double beta decay, j/g-Majorana mass, charged-current universality, e — /i — r universality, j/^ deep-inelastic 
scattering, atomic parity violation, r decays, D and K decays, Z decays, etc. Product of two couplings at 
a time have been constrained from /i — e conversion, fj, — > ej, b — > sj, B decays into two charged leptons, 
Kl ~ Ks and Bq — Bq {q ^ d, s) mass differences, etc. For a collection of all these limits see [Q. 

Let us now turn our attention to the b —>■ s"f amplitude in RPV models (Needless to add that a similar 
analysis can be carried out for b — > 1^7) Q. The aim is to generate unsuppressed diagrams for 6l — > s_r7_r. 
The RPV diagrams are shown in Fig. 1. The trilinear L- violating couplings involved are A^j-2 ^^.^ ^ij3- 
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Fig. 1: The diagrams for &i — > SB."iR induced by \[j2 3nd A^j3 couplings. 

(i) Figure la is a generic diagram in which the photon couples to the internal fermion line. It follows from 
the RPV superpotential, and as shown in the Figure la, that the fermion - scalar combinations in this case 
are (d-,L, ^il), {ujl, 1^^), and (^7^^, Mjx)- 

(ii) Figure lb is a similar to diagram where the photon is attached to the internal scalar. This time the 
fermion - scalar combinations are {viL, djh), (^i^, Uji,), and (u^l, ^Tl)- 

In Fig. 1, although chiral flips are marked only on the external b quark lines thus leading to a contri- 
bution proportional to 771^, similar diagrams can be drawn with chiral flips on the external s quark lines 
as well (leading to 6/j slJl), the latter contribution being proportional to m^. The sum of amplitudes, 
taking into account all fermion-sfermion combinations in Fig. 1, at the electroweak scale, is given by 

K'^ = - E ;"/'^^^*^ .-(m,P^ + msPn)^a,,q•'bFii , (3) 
with F]^ = QdF'^'^ + QuF'^" + Qi-F'^'^' , where 
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p^Ja ^ \-^\+ \^] , (4) 
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p.,u ^ _W_Q^ I 3^ ) + 3f G2 I 1 , (5) 
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F''' -'^G,\-t\ + 'i^G,\-^\ . (6) 

Here the Inami-Lim functions are 

Gi{x) = Ux) - Ux) .G2{x) ^ -U^) + 2Ux) - U^) , (7) 

Kn{x)-j^ l + (a;-l)z" (l-2;)"+2 ' 

The charge Q factor in the expression for F]^ reflects the charge of the internal line where the photon leg is 
attached to. The function Gi (G2) corresponds to the case that the photon attaches to the internal fermion 
(sfermion). The argument m|/m? in functions G's indicates internal lines of the fermion of generation j 
and the sfermion of generation i. Explicitly, 

2 + 5a; — a;^ a; In a; ^ , \ I — 5x — 2x^ a;2lna; 

6{i-xr +0^' ^^(") = - 6(1 -x)3 +(T^- 

Except for j = 3 in the combination A^^2-^jj3i i'^ which case there could be an internal top quark line, we 
can set the argument a: — s- to obtain 

Gi(0) = i, G2(0) = -i. (10) 

In the limit of a common mass rh for all sfermions, we then have 

2 

To be precise, we should treat the heavy top sector separately according to the detailed formulas. However, 
we use Eq. (0) only for simplicity. 




Fig. 2: The diagrams for fe^ sr^r induced by X'^gj and A^2j couplings. 

In Fig. 2, we show graphs of another type, with couplings X'^^^j and A^2j- This time we do not explicitly 
exhibit the diagrams where the chiral flips are on the b quark lines (corresponding to bj^ slJl)- They 
give rise to the amplitude 

'^i? = - E ;%^'^^^^^- .-(m,F^ + m,P,)^a,,q'^bFl^ , (12) 
Assuming degenerate masses for the left-handed and right-handed squarks, we obtain Ff^ = Q^F'^^'^. 
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Including the two types of graphs together with the SM contribution, we set up the effective Hamihonian 
responsible to the process b 57, 



w J, 



(13) 



with the chirality index h = ± (or R, L) indicating P± = ^(1 ± 75). The term involving c^h corresponds 
to the gluonic dipole contribution. At the electroweak scale the Wilson coefficients are given by the 
short-distance amplitudes outlined above, 



ij 

It is straightforward to see F^j^ = F^' and Fgj^ = F''^'^ + F'^". Also, 

ciM^w) = cf^L + E ^^3,^r2,^7k, with F^i = Q,F'^^ = F»^° 

The SM contributions are well known to be (with xt — ml/my^) 



(14) 



(15) 



7 - 5x4 - 8a;? 2xt - 3x? , 

lna:t 



I2{l-xtf 2{l-xtY 



2 + hxt — 3a;t In Xt 



4(l-a;t)3 2{l~xtY 



(16) 
(17) 



For our estimation, we evaluate the Wilson coefficients at the rub scale by the simplified renormalization 
group evolution S, 



C7h(TOfc) = ?723c7,,(tow) + I (rp3 - ry23^ C8/i(77iw) + A^f"^ , for /i = L or i?, (18) 
with 77 = as{mw) / otsiriii,) ~ 0.56. The leading log QCD corrections in the SM are given by, 



rSM 



(19) 



For the purpose of order- of-magnitude estimate, we have ignored other leading logarithmic contributions 
arising from additional operators of RPV origin, which have been outlined in Refs. p^ . 

Numerically, for the a common mass m, same for squarks and sleptons, we have 



^ ^ ij 



for E ~ 0.65 and 



D = |?7^ - 7?^ ~ -0.044 



(20) 
(21) 
(22) 



The size of the mixing-induced CP asymmetry will be controlled by the relative magnitude of Aj^ and 
Al (for the exact calculation, one should use Eq. (1)). This is given by the ratio 



R 



Ar 




C7R - 


1- {ms/mb)c7L 


Al 




C7L H 


- {ms/mb)crR 



(23) 



with the Wilson coefficients evaluated at the rub scale. 
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Upper limits 




Upper limits 


(112)(113) 


4.0 X 10-'* 


(121)(131) 


4.0 X 10--* 


(212)(213) 


2.5 X 10-3 


(221)(231) 


1.1 X 10-2 


(312)(313) 


2.5 X 10-3 


(321)(331) 


6.1 X 10-2 


(122)(123) 


8.0 X lO-"* 


(122)(132) 


3.3 X 10-3(t) 


(222) (223) 


2.7 X 10-^{*) 


(222)(232) 


3.3 X 10-3(t) 


(322)(323) 


2.7 X 10-3(=^) 


(322)(332) 


3.3 X 10-3(t) 


(132)(133) 


2.8 X lO-"* 


(123)(133) 


2.8 X 10-5 


(232)(233) 


2.5 X 10-^{*) 


(223)(233) 


3.3 X 10-3(t) 


(332) (333) 


2.5 X 10-3(*) 


(323)(333) 


3.3 X 10-3(t) 



Table 1: The existing 1-a upper limits on the A' product couplings that enter into the expressions of ctj? and cjl- The 
meaning of * and f symbols have been described in the text. A common 100 GeV mass for all scalars have been assumed 
while obtaining those limits. 

The existing upper limits (l-cr) on the magnitudes of ^'ij2^'ijz ^^'^ '^aj^aj combinations have been 
listed in Table 1 (while deriving limits the couplings have been assumed to be real in most cases) . We have 
assumed a common mass of 100 GeV for whichever scalar is exchanged. In most cases the best bounds have 
been obtained by multiplying the individual upper limits on the respective couplings, listed in Ref. |^l|. 
Note, Bg-Bg mixing (actually taking its lower limit in a conservative sense) constrains the combinations 
^'ij2''^'ij3 '^i^ sneutrino mediated one- loop box graphs ||l2|. In some cases, depending on the generation 
indices of the internal quarks and squarks, these bounds are stronger than those obtained by multiplying 
individual bounds. These have been marked (*) in Table 1. The bounds marked (f) have been derived 
from B (f)K decays and semi-leptonic B decays (Note, the authors of Ref. [|lO| have missed the color 
suppression factor in the RPV amplitudes. We have included this factor when quoting bounds in Table 1). 

As far as CP asymmetry is concerned, the size of '^o* so significant. The reason is that 

this product (even set at its upper limit) multiplied by Fy^^ adds to a much larger contribution from c^^j^ 

(see Eq. (|l5|)). The size of the other combination, namely ^j^- A^^ jA^js, is however crucial for our prediction 
for CP asymmetry. It follows from Table 1 that we can, as a rough estimate, take the magnitude of this 
combination to be ^ ±0.01. This yields rather small R ^ 3%. Such a tiny effect can hardly be interpreted 
as a signal of new physics. On the other hand, if one assumes that the squarks are much heavier than 
the sleptons, as happens in the gauge- mediated supersymmetry breaking models ||l3|| , then in the RPV 
amplitudes one can consider only the slepton- mediated graphs, as the squark- mediated ones sufficiently 
decouple to leave no numerical impact. For our remaining discussions, we stick to the latter scenario. In 
this case the factor of ^ in Eq. ( ^l| ) becomes 1 and the variable D in Eq. (|2^) is replaced by 

D' = ^r;M _ |,yM ~ 0.25 , (24) 

which is rather large compared to D. This time, with | ^ Kj2'^ij3 \ ^ 0.01, we obtain R as large as 16%. 
Note that if we set this combination to be 0.04 (0.05), which is not very unrealistic, as the limits have been 
set only on an order-of-magnitude basis, R may go up to ~ 70% (90%), which is a rather large effect. 

Although in the case of CP asymmetry the impact of \2j^'aj is insignificant as noted earlier, this 
combination, both in magnitude and phase, may turn out to be quite significant in the total branching ratio 
(BR) for h S7. Since the BR is proportional to (cy^ -I- c^^), the RPV contribution to cjl in Eq. (|2l| ) 
interferes with a rather large SM contribution and hence appears linearly, while the RPV contribution to 
C7R (being the only contribution) appears only quadratically. Also, within cyl, the interference between the 
SM and the RPV pieces may be either constructive or destructive depending on their relative phase. The 
possibility of this partial cancellation may enable us to admit larger values of RPV couplings required to 
generate a sizable CP asymmetry, while stay perfectly consistent with the BR constraints. In Fig. 3, as an 
illustrative example, we exhibit our prediction for the BR as a function of jA^j, where A^ = ■ \ij2\j3 = 
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KsjKlj- The different angles refer to the relative phase between the two pieces in Eq. In Fig. 4, 

we demonstrate how large CP asymmetry (only its absolute value) can be generated as a function of |A^| 
without violating the BR constraint. For the latter, we used Eq. (1) with the sine of the phase combination 
set to unity. Here we wish to make the following remark. While drawing the new physics curves in Fig. 3, 
we took a conservative approach that the SM reference line is at its present mean value. One might as well 
place the SM reference line close to its lower limit so that larger values of |A^| could be allowed leading to 
a possibility of larger CP asymmetry. 

Before closing, we highlight the following salient features: 

1. The sign of the RPV induced CP asymmetry is arbitrary, since the signs of the A' couplings are a 
priori unknown. It is therefore possible to have a CP asymmetry not only large compared to the SM 
prediction but also with an opposite sign. The occurance of a sign flip, in particular, may constitute 
an unmistakable signature of new physics . 

2. In principle, the A' couplings need not be complex in order to generate CP asymmetry. More precisely, 
in Eq. (1) even with (pn = (i.e. with real A' couplings), one may obtain a non-zero sine function 
from the B-B mixing phase (f)M or the SM decay phase (j)^. Nevertheless, a non-zero 0/? may trigger 
a much larger CP violating effect. 

3. Here we present a simple way of understanding how large CP asymmetry can one predict without 
violating the branching ratio constraint. The present experimental BR on B Xgj is (3. 15 ±0.54) x 
10""* Ijl^l, while the present SM next-to-leading order estimate for BR (6 57) is (3.28 ±0.33) x 10^^ 

To check whether our prediction of a large CP asymmetry at all contradicts the constraints 
on the BR, let us consider the following illustrative example. Assume that new physics contributes 
only to ctr. In our case, this amounts to saying that Kj2Kj3 non-zero, while K3j'^i2j 
is vanishingly small. Recall, to a good approximation, cth/c-jl = a is a quantity that controls the 
size of the CP asymmetry (for an exact expression, see Eq. (1)). Then the BR, being proportional 
to {c^ji + c^l), is modified by an overall factor (1 + a^). Now to demonstrate how large a one can 
tolerate, let us imagine a situation when the SM prediction of the BR is at its present 90% CL lower 
hmit, i.e. 2.74 x lO^'*, and the experimental BR settles at its 90% CL upper limit, i.e. 4.04 x 10^*. 
This means that within the current experimental and theoretical uncertainties a 50% (= a^) new 
physics induced enhancement of BR is perfectly tolerable, and hence a CP asymmetry as large as 
2a/ (1 + a^) ~ 100% (putting a = 0.7 in Eq. (1)) is very much consistent with the present constraints 
on the BR at 90% CL. 

4. The RPV models are better placed than the left-right symmetric model in generating large CP asym- 
metry. The reason is that the R-parity conserving MSSM contribution to cyl may interfere de- 
structively with the SM amplitude in some region of the parameter space [|l^. Just like the RPV 
contribution to ctl , this also helps to keep the BR under control while admitting large RPV couplings 
responsible for large CP asymmetry. More specifically, the charged Higgs (chargino) loops interfere 
constructively (destructively) with the SM W loops, and we know that in the exact supersymmetric 
limit the net amplitude vanishes p8[ |. 

5. The combinations ^ii2\i3, within their current experimental bounds, not only induce large CP asym- 
metry in 6 ~> S7 decay, but also enhance the direct CP asymmetry in B^ tt^K channels jl^. Hence 
if upcoming B factories indicate large CP asymmetry in both these (uncorrelated) channels, and per- 
haps with signs opposite to their respective SM predictions, it may constitute an unambiguous signal 
of the RPV scenario. We must admit though that the extraction of 7 from B ttK would be even 
more difficult in the presence of new physics contamination. 

GB thanks the following centers for hospitality where portions of the work were done: NCTS, Hsinchu, 
Taiwan; CERN Theory Division, Geneva, Switzerland; and MPI Kernphysik, Heidelberg, Germany. DC 
wishes to thank SLAG Theory Group and U.S. DOE for hospitality while revising the final manuscript. 
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Fig. 3: Br(fe —^ sy) versus |A^| for various choices of the relative phase between KajKij and VaVf^. The SM 
prediction and the experimentally allowed region are shaded. 
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Fig. 4: The magnitude of the CP asymmetry for various choices of the relative phase between X^ij^W^iSj Vt(,V(*. 
Here we assume that both the sine factors are unity in Eq. (1). 



9 



